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Abstract. 

Photoionization cross section of the fullerene ion Cqq has been calculated within 
a single-electron approximation and also by using a consistent many-body theory 
accounting for many-electron correlations. Electronic wave functions of the ground 
and excited states have been determined within the jellium model based on the local 
density approximation. It is predicted a giant resonance in the photoionization cross 
section at the photon energy ^ 25 eV. It is demonstrated that the resonance profile is 
much higher and narrower than in the case of the neutral fullerene Cgo- 

Since the discovery till nowadays the fullerene Cqq attracts a permanent rapt 
attention of investigators. From the physical point of view an interest to this object 
is associated with its exotic hollow, but highly stable structure caused by a significant 
delocalization of a large number of valence electrons in this molecule. These features 
manifest themselves in a number of interesting physical phenomena occuring in collision 
processes involving fuUerenes. For example, surface plasmon oscillations in the fullerene 
Ceo can be excited. The plasmon excitations have a profound collective nature and 
they influence significantly the formation of the cross sections in collisions of various 
kind involving fullerenes p]. The plasmon excitations in the fullerene Ceo have been 
well investigated both experimentally and theoretically |21ISllllElinilZlElinilIIIllIIllI2!. 
Usually, they manifest themselves as giant resonances in the excitation spectrum of Ceo- 
In this connection, there arises a question on the possibility of excitation of the plasmon 
oscillations in fullerene Ceo ions of different multiplicity. Recently, the photoionization 
cross section of singly charged fullerene ion C^ has been measured ^31- However, 
theoretical calculation of this cross section has not been performed so far. 

In this work, for the first time, we have calculated the photoionization cross section 
of the fullerene ion C^q within the photon energy range from the Cqq ionization threshold 
up to 40 eV. This energy range has been chosen, because the photoionization cross 
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section for the neutral fullerene Ceo possesses a very strong giant resonance centered 
at the photon energy about 20 eV [HI ■ Using the theoretical framework developed 
earlier by us and applied for the description of photoionzation of the neutral fullerene Ceo 
we have determined the location of the giant resonance in the photoionization 
spectrum of Cqq (to be equal ~ 25 eV) and described it's shape. 

The atomic system of units is used throughout the paper, me = \e\ = h = 1. 

The total photoionization cross section is equal to a sum of partial photoionization 
cross sections for each of the fullerene ion orbitals. The amplitude of ionization of 
each orbital is calculated both within the single-electron approximation and also by 
taking into account many-electron correlations. At the first stage, we use the frozen 
core model and the local density approximation assuming that there is a single electron 
transition only during the photoionization process. Then, the correlations between the 
transitions from different electron states are taken into account within the random phase 
approximation (RPA). The RPA amplitudes of photoionization Di^^uii^) are found from 
the matrix form of the RPA equations flM^ 

U3U4 

Here, the matrix elements D^^ui are expressed via the single-electron photoionization 
amplitudes du2Ui, the matrix elements of virtual electron-hole excitation xi^)^ and the 
matrix elements of interelectron interaction Vu^i^j^i,^,^^ Note that the excitation 

spectrum has the Rydberg series of discrete excited states because the self-consistent 
potential of the fullerene ion has the Coulomb asymptotic. Therefore, the virtual 
electron transitions to the discrete excited states of the fullerene ion play more important 
role than in the neutral fullerene Ceo case. 

The photoionization cross section and the oscillator strengths calculated both 
within single-electron approximation and also by taking into account many-electron 
corellations have been checked on the consistancy with the sum rule. 

The inter-electron interaction within the local density approximation reads as 
follows [12] 

1 d^etTln] 



6{f-f'), (2) 

n=no (r) 



|r*— r*'! dn 

where e^°'"[n] is the exchange-correlations energy of homogeneous electron gas with 
density n within the Gunnarson-Lundqvist approximation ^H] 

3/9^/' 1 



, ^ , , - 0.0333G r, (3) 

G(x) = (l + x3)ln(l + i)-x2 + ^-l, (4) 

1 /3 

where rs{r) = (^ ^^^^^^ ^ is the Wigner-Zeith radius, no is the equilibrium density. 

Single-electron states of the fullerene ion have been calculated within the spherical 
jellium model with the self- consistent potential fTH. IT^ determined within the local 
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density approximation. Analogously to an atom, the single-electron wave functions of 
Cqq within this model, are given by 



where rii, k, rrii is the usual set of quantum numbers for the i-th spherical shell; Yim{0, (p) 
is the spherical harmonic; Pn^i- (r) is the radial part of the wave function determined in 
the self-consistent potential of the fuUerene ion. The number of nodes of the radial wave 
function P„j-(r) is equal to (nj — — 1). The main limitation of this method arises from 
the electrons self-interaction. Note that the self-interaction correction does not change 
the order of the electronic energy levels and the value of the total electronic energy, but 
it alters the value of the ionization potential very significantly. 

As in the case of neutral fullerene, we assumed that four 2s2p electrons of each 
carbon atom are delocalized in the fullerene ion Cqq. The ionic core of the fullerene is 
formed by fourfold-charged positive carbon ions C^"*" with the electronic configuration 
Is^. The charges of the carbon ions are averaged over the sphere of radius R, being 
a hypothetic radius of the fullerene ion. In this paper we neglect the variation of the 
fullerene's radius caused by its ionization and put the radius of the fullerene ion equal 
to the radius of the neutral fullerene, i.e. R = 6.69 a. [TB] . 

For better understanding of the electronic properties and interatomic bonding in 
the fullerene molecules, it is useful to establish the correspondence of the electronic 
eigenstates introduced in our model to the a and vr orbitals widely used for modeling 
the planar graphite surface ^Tj. The ground state configuration for the ion C^q has been 
chosen in a similar way as for the neutral fullerene Cqq jTT] . Removing one electron from 
the highest occupied molecular orbital of Ceo, one derives 



Here an atomic notation of electronic configuration n is used . 

The calculated energy level spectrum is shown in figure 1. The ionization potential 
5.5 eV obtained within the LDA differs from the experimental value 11.5 eV ^H]- The 
origin for this difference is discussed above. 

The maximum number of electrons in a shell characterized by the orbital momentum 
/ is equal to Ni = 2{2l + 1). The most outer a— (1=9) and vr— shells (1=5) in Cqq are 
open and have A^e = 18 and A^e = 9 electrons respectively, whereas A^g = 38 {a— shell) 
and = 22 (vr— shell). This means that the total orbital moment L and the total spin 
5* of such a system are non-zero. 

It is well known that an accurate calculation of the photoionization cross section of 
an atomic system with open electronic shells is much more complicated than that for a 
closed shell atomic system, because for an open shell system the calculation involves a 
huge number of additional terms which vanish in the case of a closed shell system |19j . 
In order to avoid this technical problem we have used the averaged term approximation. 



0i = 0n.,ii,m,(r) = -Pn,lXr)Yi^mA^,'^) 



(5) 
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Figure 1. Calculated energy levels spectrum of the fuUerene ion Cqq. 
none-nodes and one-nodes orbitals are partially filled. 



The outer 



We have assumed that the electrons from the most outer open shells behave in the 
photoionization process as if the shells would be closed. This implies that one can put 
L = 0, 5 = when calculating the partial contribution to the photoionization cross 
section from these shells and take into accout their filling ratios gi = Ne/2{21 + 1) < 1. 

It is natural to expect that the averaged term approximation affects the 
photoionization cross section behavior. The numerical analysis shows that the cross 
section calculated within the single-electron averaged term approximation does not obey 
the sum rule. Therefore, one should introduce appropreate correction when calculating 
the photoionization cross section in the random phase approximation. Thus, we have 
input extra coefficients decreasing the amplitude of electron-hole excitations from the 
open shells. These coefficients have been determined from the sum rule analysis of 
the partial photoionization cross sections of the open fuUerene shells performed within 
the single-electron approximation. The coefficients enter as numerical factors in the 
matrix elements of operator x{^) equation (0). Although this procedure improves 
the consitence with the sum rule of the photoionization cross sections calculated within 
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Figure 2. Photoionization cross section of tlie fullerene ion Cqq calculated within 
different approximations. Dotted line is the result of RPA and the averaged term 
approximation, dashed line is the same as dotted line but obtained with the optimal 
values of additional coefficients in the averaged term approximation, solid line is the 
same as dashed line but normalized according to the sum rule. 

the RPA, it leaves the deviation of the RPA photoionization cross sections from the sum 
rule on the level of 30 %. 

Figure 2 shows the Cqq photoionization cross sections calculated by different 
methods: RPA and the averaged term approximation, RPA with extra coefficients, RPA 
with extra coefficients normalized on the sum rule. The last curve is the most accurate 
result of our calculations of the Cqq photoionization cross section. These calculations 
demonstrate that the giant resonance pattern is not affected much by the choice of the 
coefficient in equation dT}, but the maximum value of the cross section depends strongly 
on their choice. 

The calculated cross section is compared with theoretical pT^ and experimental 
cross sections jH] for the neutral Cqq in figure 3. In order to elucidate the difference 
between the photoionization cross sections for Ceo and Cqq, the Cqq photoionization 
cross section is shifted in figure 3 towards the ionization threshold of Cgo- Under the 
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Figure 3. Comparison of tlie photoionization cross sections for the neutral fullerene 
Cgo and the ion C^q. Solid line is the most accurate result of calculation of the 
photoionization cross section for the fullerene ion C^q. This curve is shifted on 

A = 7^60 - /p«", where Ip's are the ionization potentials. Dashed line shows the 
experimental photoionization cross section for Cgo IHl j dotted line is the photoionization 
cross section for the neutral Cgo calculated within our model [12]. 



condition of coincidence of the ionization thresholds of Ceo and C^q, the giant resonance 
in the Cqq photoionization spectrum is placed at the photon energy 21 eV. It is well 
known that the resonance energy of the plasmon oscillation can be easily evaluateded 
for a spherical fullerene with the use of the classical Mie formula 



\ {21 + l)/?3 ' W 

where uji is the frequency of a surface plasmon oscillation with angular momentum / and 
Nf. is the number of delocalized electrons. For the diploic plasmon mode, / = 1, in the 
fullerene Cqq with R = 6.69 a.u. [12] and = 240, one derives Ui = 20 eV. 

Since we have assumed the same radius for the fullerene ion as for the neutral 
fullerene and there is only a small difference in the number of delocalized electrons, it 
is natural to expect that the plasmon resonance energies in the ion and in the neutral 
fullerene are close. Figure 3 illustrates this fact. It shows that nevertheless there is a 
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slight shift of the position of the giant resonance in the photoionization spectrum, 
which might be a result of a slight variation of the radius of the fuUerene ion neglected 
in our calculations. 

Our analysis demonstrates that the giant resonance in the photoionization cross 
section for the fuUerene ion arises due to the correlation between transitions from a— 
orbitals. This fact proves the similarity of the nature of the giant resonances in the 
cases of the neutral fuUerene and the ion. Qualitatively, one can explain the origin of 
the giant plasmon resonance by oscillations of the a— electrons density against the ionic 
core. 

It is also necessary to point out that the giant plasmon resonance in the 
photoionization spectrum of Cgg is much stronger as compared with that for Cgo- The 
similar behaviour is well known from the investigation of giant resonances in many- 
electron atoms and their ions. It can be explained by enhanced influence of the Coulomb 
core potential on the ground and excited states of valence electrons in the case of ionic 
target 

In conclusion we stress that the method developed in our work allowed us to predict 
the plasmon resonance pattern in the photoionization cross section of the fullerene 
ion Cqq. This method can be straightforwardly utilized for the calculation of the 
photoionization cross sections of positive and negative ions of Cqq of larger multiplicity, 
which are of great interest for experimental investigations. 

The authors acknowledge support of this work by the INTAS (grant No 03-51- 
6170) and the Russian Foundation for Basic Research (grant No 03-02-16415-a). 
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